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We analyze theoretically the oscillations that the magnetoresistivity of two- 
dimensional electron systems present when a high intensity direct current is applied. 
In the model presented here we suggest that a plasma wave is excited in the system 
producing an oscillating motion of the whole two-dimensional electron gas at the 
plasma frequency. This scenario affects dramatically the way that electrons interact 
with scatterers giving rise to oscillations in the longitudinal resistivity. With this 
theoretical model experimental results are well reproduced and explained. 
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Magnetotransport properties of highly mobile two-dimensional electron systems 
(2DES) is a subject of increasing interest. The future nanoelectronics will integrate 2DES 
as key parts of different devices as transistors, diodes or capacitors. In particular the phe- 
nomena of microwave-induced resistivity oscillations (MIRO) and zero resistance states 
(ZRS^jMiSiSiSi have attracted considerable attention both from theoretical and experi- 
mental sides. However the mechanism responsible is still under debate. 

Another striking non-linear effect has been observed when a 2DES is subjected to a 
high intensity DC-current&i&ii. These experiments report oscillations in the longitudi- 
nal magnetoresistivity (p xx ) that are periodic in inverse magnetic field and tend to be 
qualitatively similar to the ones obtained in the experiments of microwave (MW) excited 
magnetoresitivity in Hall bars^ 1 ^. Another remarkable experimental outcome is that 
under DC-current excitation, the p xx response shifts towards higher magnetic fields (B) 
with increasing DC-current intensity. Also the electron sheet density dependence shows 
that the maxima positions in the p xx oscillations are scaled as 1/ ^/nl where n e is the elec- 
tron sheet density. The first explanation to these effects*^ 4 ^ considers that oscillations 
are related to elastic scattering mediated by a short range disorder between Hall field- 
tilted Landau levels. Accordingly an electron may transfer momentum in the x-coordinate 
that is equivalent to a hopping in the y direction. p xx reflects this hopping through the 
conductivity a yy . Other theoretical proposals have been reported recently^. 

In this letter we present an alternative approach to explain these effects. In our model 
the DC-electric field (Edc), which produces the DC-current (Idc), displaces the orbit 
center of the harmonic quantum oscillators that are the electrons in the presence of a 
uniform and perpendicular B. This displacement with respect to the positive lattice ions 
disturbs the 2D electron gas leading to a collective excitation of the system. In other 
words, a plasma wave originates from a spatial shift of the 2DES as the charge seeks to 
move to restore equilibrium. The system responds with an oscillatory motion where the 
frequency is denoted as w p or plasma frequency. Evidently its effect is only measurable 
at sufficient high DC-current intensities. This scenario affects dramatically the way that 
electrons interact with scatterers and is reflected in the p xx response. We explain and 
reproduce the p xx oscillations periodicity in B~ l , the shift to larger B with increasing 
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current intensity and the dependence on n e . 

The electrons of a 2DES subjected to a perpendicular B, behave as harmonic quantum 
oscillators. Then if we apply a constant electric field Ejjc i n the current direction (x- 
direction), we obtain a displacement x of the center of all oscillators in the same direction: 
the entire 2D electron gas moves. The maximum value of x is given by^>^ x max = 
where e is the electron charge and w c the cyclotron frequency. The spatial shift of the 
electron gas as a whole with respect to the fixed positive background (lattice ions) give rise 
to two lines of opposite charge at either end of the sample creating an electric field E p (see 
Fig. 1). We obtain this field using basic Electrostatics 15 : E(p) = 2 2tT f * being L e f the 
average effective distance between the two lines of opposite charge that we approximate 
by L e f = L — Xmax-, e is the dielectric constant and L the sample length. This field tends 
to restore the system to its equilibrium position producing in the system a collective 
excitation or plasma wave. Eventually the 2D electron gas obey the equation of motion 
of a harmonic oscillator 16 that leads to oscillations at the plasma frequency w p : 



On the other hand the plasma oscillation can be damped. Physically the coupling of 
the plasma to the lattice through phonon scattering results in a net decrease of the 
energy of the oscillating electron gas. The result is a progressive decay of the plasma 
oscillation being reflected in a continuous decrease of its amplitude. Eventually the plasma 
oscillations can be totally damped. However the external field E^c is still acting on the 
2DES exciting the plasma wave and avoiding the oscillatory motion to stop. Thus Ejjc 
plays the role of and exciting force on the oscillating plasma. If we name this force as 
Fjjc, the plasma equation of motion can be written as: 

Nm *~T + Nm*>y^- + Nm*w 2 p x = F DC (2) 

where N is the number of electrons and 7 is a phenomenologically introduced damping 
factor which accounts for the electronic interactions with the lattice yielding acoustic 
phonons. To solve this differential equation we propose as a possible solution a harmonic 
function like x p = Aoc cos w p t. The expression of Fdc an d the amplitude Aoc are to be 
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determined. Thus we obtain: 

eE DC ( , 

x v = cosw p t Id) 

if Fdc = NeEf)c cos wt and Adc = ^*°^ ■ Therefore we have a 2DES that oscillates 
as a whole with w p , displacing the electrons orbit center with x p (t). The oscillations are 
excited by Fdc an d damped by coupling with the lattice. 

Now we introduce the scattering suffered by electrons due to charged impurities ran- 
domly distributed in the sample. Without high intensity DC-current excitation, an elec- 
tron in an initial state corresponding to an orbit center position scatters and jumps 
to a final state with orbit center X^, changing its average coordinate in the static electric 
field direction in AX° = X^ — X®. This is the relevant direction (x-direction) in this 
problem due to the presence of the DC electric field. It means that all the scattering 
jumps will be averaged to zero except the ones in the x-direction. Under a high intensity 
DC-current the plasma wave is excited and the electronic orbit center coordinates change 
being given, according to our model^ 1 ^, by: 

X p = X° + x p = X° + A DC cos Wpt. Thus, due to the collective oscillation all the electronic 
orbit centers in the sample oscillate back and forth in the x direction through x p . When 
an electron suffers a scattering process with a probability rate given by W m>n , it takes a 
time r = — for that electron to jump from an orbit to another. This probability rate 
is calculated according to the model described in ref^ 4 -. If we consider that the electron 
jumps from the initial state oscillation middle position, and it takes a time r to get to 
the final one, then we can write for the average coordinate change in the x direction: 
AX P = AX° + Adc cos w p T Finally the longitudinal conductivity be calculated: 

&xx oc J dE^—(fi — ff), being /j and // the corresponding electron distribution func- 
tions for the initial and final states respectively and E energy. To obtain p xx we use the 
well-known tensor relation p xx = 2 a T 2 1 where a xy ~ 

In Fig. 2a, we present calculated p xx versus B at a direct current Idc — 300/iA We 
observe clear oscillations in p xx showing four peaks. The oscillations can be explained in 
similar terms as the ones obtained with MW excitation^(see Fig. 3). When no plasma 
wave is excited, electrons jump between fixed orbits and on average an electrons advances 
a distance AX° (Fig. 3. a). When the DC-current excitation is on, orbits are not fixed 
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and instead move back and forth through x p with the frequency w p . When the orbits, 
due to plasma oscillation, are moving backwards during the scattering jump the electrons 
advance an average larger distance than the no plasma wave case: AX P > AX°. This 
corresponds to an increasing conductivity (Fig. 3b). When the entire electron gas is 
moving forward, during the jump the electron advances on average a shorter distance: 
AX P < AX° (Fig. 3c). This corresponds to a decrease in the conductivity with respect 
to the case without DC-current excitation. 

In Fig. 2b we present calculated p xx versus We observe that p xx is periodic in 

B^ 1 in agreement with experiments^^ with an spatial period of 5, (see figure). This 
is clearly shown in the inset where we present the inverse of peaks position, B~ x versus 
the order of the peaks n. According to our model, p xx oc cosw p t = cos^ 13 , K being a 
constant, i.e., p xx is periodic with B" 1 . Thus, in the peaks the next condition is fulfilled: 
B^ 1 = ^-2im. Therefore we obtain the equation of a straight line that crosses the origin 
as in experiments^ 1 ^. 

In Fig. 4, we present calculated p xx versus B for different values of Ibc- We observe 
that the peaks shift towards higher B with increasing Idc- An increasing Idc corresponds 
to a larger Edc an d x max affecting eventually w p (see eq. 2) that becomes larger too. 
Considering the functional dependence of p xx through cos w p t, a larger w p give rise to 
more peaks and a shift of p xx response to larger B. A similar behavior is obtained in 
MW-excited p xx response in Hall barsi. In Fig. 5, we present the calculated peaks 

1/2 

position B n versus n e . Two straight lines crossing the origin are obtained for peak 
order 1 and 2. The inset shows p xx versus B for different n e . Again following our model, 

1 /2 

p xx oc cos rJf p = cosC^^-b = cos ,,~ — , C being a constant. Peaks maxima fulfill 
2im = y 2 =^> B = -^{-Tj2- As in the experimental outcome^, we obtain the equation 
of a straight line crossing the origin. 
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Fig.l caption: Schematic diagram showing electronic transport through a Hall Bar. 
A high intensity DC-current produces a displacement x max of electrons orbit center with 
respect to the fixed positive background of lattice ions. The resulting linear distribu- 
tions of opposite charge in opposite sides of the 2D sample give rise to an electric field Ep. 

Fig.2 caption: (a) Calculated p xx versus B at a direct current Idc = 300/iA We 
observe clear oscillations in p xx response. Peaks order from 1 to 4 is shown, (b) 
Calculated p xx versus inverse magnetic field (-B -1 ). The inset shows that p xx is periodic 
in being 5 the spatial period. T=1K. 

Fig. 3 caption: Schematic diagrams of electronic transport through the 2D sample 
without (a), (fixed orbits) and with plasmonic excitation (b) and (c), (oscillating orbits). 

Fig. 4 caption: Calculated p xx versus B for different values of I DC (50pA — > AOOpA). 
We observe, as in experiments, that the peaks shift towards higher B with increasing 
Idc- The p xx shift and the increasing number of peaks with Idc suggest and increasing 
w p . T=1K. 

Fig. 5 caption: Calculated peak position B n versus nj 2 . n e is the electron sheet 
density. Two straight lines crossing the origin are obtained for peak order 1 and 2. The 
inset shows p xx versus B for different n e . T=1K. 



